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f*^ i We revise the stability of the tracking solutions and briefly review the potentials of 

VQ , quintessence models. We discuss the evolution of linear perturbations for V(<j>) = 

^h ■ Vb exp(A</< 2 /2) potential in which the scalar field is non-minimally coupled to cold 

^S ' dark matter. We consider the effects of this coupling on both cosmic microwave 

VQ | background temperature anisotropies and matter perturbations. We find that the 

f *^ , phenomenology of this model is consistent with current observations up to the 

"""■^ ■ coupling power n c < 0.01 while adopting the current parameters measured by 

i-pH ■ WMAP, £^ 0) = 0.76, ClfJ m = 0.191, Q { b 0) = 0.049, and h = 0.70. Upcoming 

cosmic microwave background observations continuing to focus on resolving the 

/"} ■ higher peaks may put strong constraints on the strength of the coupling. 
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1. Introduction 

If we treat Type la supernovae (SNe la) as standardized candles, then 
the Hubble diagram of them shows that the expansion of the Universe is 
currently accelerating 1 . Combining measurements of the acoustic peaks in 
the angular power spectrum of the cosmic microwave background (CMB) 
anisotropy 2 and the matter power spectrum of large scale structure (LSS) 
which is inferred from galaxy redshift surveys like the Sloan Digital Sky 
Survey (SDSS) 3 and the 2-degree Field Galaxy Redshift Survey (2dFGRS) 
4 has also confirmed that a component with negative pressure (dark energy) 
should be added to the matter component to make up the critical density 
today. 

A quintessence field is a dynamical scalar field leading to a time de- 
pendent equation of state (EOS), u)$. The possibility that a scalar field at 
early cosmological times follows an attractor-type solution and tracks the 
evolution of the visible matter-energy density has been explored 5 . This 
may help alleviate the severe fine-tuning associated with the cosmological 
constant problem. However we need to investigate the tracking condition 
and its stability at the matter dominated epoch carefully 6 . 
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Are there experimental ways of checking for the existence or absence 
of dark energy in the form of quintessence? There are several different 
observational effects of matter coupling to the scalar field on CMB spectra 
and matter power spectrum compared to the minimally coupled models . 
And this can be used to check the existence of quintessence. 

This paper is organized as follows. We briefly investigate the con- 
dition and the stability of tracker solution and review the potentials of 
quintessence models in the next section. In Sec. Ill, we show the cou- 
pling effects on CMB and matter power spectrum. We conclude in the last 
section. 

2. Quintessence Models and Tacker Solutions 

Many models of quintessence have a tracker behavior, which solves the 
"coincidence problem" (i.e. initial condition). In these models, the 
quintessence field has a density which closely tracks (but is less than) the ra- 
diation density until matter-radiation equality which triggers quintessence 
to start behaving as dark energy, eventually dominating the Universe. This 
naturally sets the low scale of the dark energy. However the present small 
value of dark energy density still cannot be solved with quintessence ( "fine- 
tuning problem"). 

Since the energy density of the scalar field generally decreases more 
slowly than the matter energy density it appears that the ratio of the two 
densities must be set to a special, infinitesimal value in the early Universe in 
order to have the two densities nearly coincide today. To avoid this initial 
conditions problem we focus on tracker fields. 

2.1. Tracking Condition 

We rely details of this section on the reference 6 due to the shortage of 
space. We introduce new quantity 9 which is related to the ratio of the 
kinetic energy and the potential energy of the scalar field : 

2tf = ln^=lni±^ (1) 

PE l-w w 

9 can have any value and especially positive 9 means the kinetic energy 
dominated era and negative 9 indicates the potential energy dominated 
one. Now we can define EOS as the function of this new quantity 9 : 

u)$ = tanh# (2) 
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The "tracker equation" (T = V"V/(V') 2 ) can be expressed by 9 ; 

_ 3 {Auj r - cj ) (1 - n^) 1-UJ0 1 9 

2(1+ W4,) (3 + 9) 2(1 + wt) (3 + 9) (1 + c^ ) (3 + 9f U 

where uv is EOS of the radiation, A = a eq /(a + a eq ), a eq is the scale factor 
when the energy density of the radiation and that of the matter become 
equal, and tilde means the derivative with respect to x = In a = — ln(l + z) 
a . This equation looks like quite different from the well known tracker 
equation. But when we choose the early Universe constraints (A ~ 1 and 
n^ ~ 0) we can get the well known tracker equation. 
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When we can ignore the change of 9 (i.e. when lo^ is almost constant), we 
can get the tracking solution. From the above tracking equation we can 
check this condition. 

r* i + i^^(i- n,) (5) 

l 1 + U)<p 

This equation can be rearranged to see the behavior of the EOS as following. 

u r A(i-n <t> )-2{r-i) 
^ " (i-^) + 2(r-i) (6) 

To investigate this equation more carefully, we define the new quantities. 

Q = (1 - fy) (7) 

f = (r - 1) (8) 

where Q shows the energy information of the Universe and F depends on 
the form of the given potential. With these we can rewrite the equation 
(6). 

Q * 2F 

loa, ~ qv— A — — (9 

Q + 2F Q + 2F y ' 

In the reference 5 , this equation is expressed as : 

u> r -2F 
Ud> ^ — (10) 



a Where we put the present value of scale factor, a*- ' 
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This equation (10) can be true only when Q, A ~ 1, which can be satisfied 
at the early Universe and not at the late one. So with this equation, it 
is not proper to check the evolution of tracking solutions at late Universe. 
Instead we should use the equation (9) to check the evolution of the tracking 
solutions. Before checking the properties of this equation, we should notice 
that Q is always positive and has the interval as < Q < 1. F can be 
positive or negative based on the given shape of the potential. 



2.2. Stability Of Tracker Solution 

We need to check that solutions with w<£, which is not equal to the tracker 
solution value (wo) can converge to the tracker ones (i.e. Are tracker solu- 
tions stable?). To check this we need to check the small deviation (Suj) of 
the tracker solution of EOS. 

cj0 = uj + Slj (11) 

If we insert this into the tracker equation (3), then we have following. 



5u 



3r 
"2 



9. 



{Au> r -u>o)(l-n<f,)+(l-uJo) 6u+-(1-u>q) {l+Au> r ){l-{l<f,) 



Slo = 
(12) 

where we use the tracking condition (5). The general solution to this nonlin- 
ear differential equation cannot be obtained analytically. But this equation 
can be simplified as follow in the early Universe constraints. 



~ 3r 



(1 + w r ) — 2ujq 
The solution of this equation is 



Suj + — (1 — u>q)(1 + Lu r )5tu ■ 







Suj 



,7i 



(13) 



(14) 



where 



7i 



3rl 



(1 + w r ) - Wo 



± 



rl 



18(1 + w r )(l - LU ) - 9 -(1 + Wr) - ^0 



(15) 

The real part of this is negative for oj less than 2/3. So Su> will decays 
exponentially and solution reaches to the tracking one. In addition to this 
it also oscillates due to the second term. For the late Universe case, we can 
change this equation. In the late Universe A goes to zero and Q,p is not 
zero. With these we can modify the general equation (12). 



~ 3 r i~9 

Suj +- (1 + f^Wo) - 2w Suj + -(1 - Wq)(1 - il^Suj ■ 







(16) 
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Table 1. Quintessence models. 



Potential 



Reference 



Properties 



Vb exp (-A$) 

V /(j> a ,a>0 

m 2 <t> 2 ,\(j> 4 

Vb(exp M p /4>- 1) 

Vbcxp(A0 2 )/<r 

Vb(coshA</>- 1)p 

Vb sinh-" (Xcp) 






Vb(e 
Vb[(</>-B) 

Vb expp^/Mj,) 2 
Vb cosh[A<£/Mp] 



Ratra & Peebles (1988), Wetterich (1988) 5 

Ferreira & Joyce (1998) 5 

Ratra & Peebles (1988) 5 

Frieman et al (1995) 8 

Zlatcv, Wang & Stcinhardt (1999) 9 

Brax & Martin (1999,2000) 10 

Salmi & Wang (2000) n 
Sahni & Starobinsky (2000) 12 , 
Ureiia-Lopcz & Matos (2000) 13 
Barreiro, Copeland & Nuncs (2000) 14 
Albrccht & Skordis (2000) 15 
Lee, Olive, & Pospclov (2004) 16 
Lee, Olive, & Pospclov (2004) 16 



w = A 2 /3- 1 

A > 5.5- 4.5, il < 0.1-0.15 

o> > -0.7 

PNGB, M 4 [cos(0//) + 1] 

f7 m > 0.2, u) < -0.8 

a > ll,w ~ -0.82 

p < 1/2, o) < -1/3 

early time : inverse power 

late time : exponential 
a > 5.5,/3 < 0.8,0) < -0.8 

O) 1 

O) 1 

o>~-l 



We can repeat the similar step to find the solution of this equation if we 
assume that Qj, is almost constant. 



5lu oc a 12 



(17) 



where 



72 



1 1 z / rl I 2 

-(1+£Vjo)-woJ±2 y 18(1 + 0/0(1 - wo) - 9[-(l + O^o) - w 



(18) 



The real part of this solution can be negative if wq satisfies following. 



w < 



(2-fy) 



(19) 



where 1 < (2 — £1^) < 2 for the entire history of Universe. 



2.3. Quintessence Potentials 

We display the potentials of the quintessence models in Table 1. Any detail 
of each model can be found in each reference. 
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3. Coupled Quintessence 

The general equation for the interaction of a light scalar field tp with matter 
is, 

S, = /dW^j^Vw -R}- V(<f>) BzMfWfWv 

+ YfiiWi - B J WmjWj] \ • (20) 

The coupling gives rise to the additional mass and source terms of the evo- 
lution equations for CDM and scalar field perturbations. This also affects 
the perturbation of radiation indirectly through the background bulk and 
the metric perturbations 7 , 17 . The value of the energy density contrast of 
the CDM (f2 c ) is increased in the past when the coupling is increased. We 
specify the potential and the coupling as in the reference 7 , 16 . 

VM = V u ex P (^), exp[B c (*)] = i ^^ Y (21) 






3.1. CMB 



Now, we investigate the effects of non-minimal coupling of a scalar field to 
the CDM on the CMB power spectrum. Firstly, the Newtonian potential 
at late times changes more rapidly as the coupling increases. This leads to 
an enhanced ISW effect. Thus we have a relatively larger Cg at large scales 
(i.e. small £). Thus, if the CMB power spectrum normalized by COBE, 
then we will have smaller quadrupole 18 . This is shown in the first panel 
of Figure 1. One thing that should be emphasized is that we use different 
parameters for the ACDM and the coupled quintessence models to match 
the amplitude of the first CMB anisotropy peak. The parameter used for the 
quintessence model is indicated in Figure 1 (i.e. tt[ = 0.76, flm = 0.191, 

Q b = 0.049, and h — 0.7, where h is the present Hubble parameter in the 
unit of 100kms _1 Mpc _1 ). However, these parameters are well inside the 1 a 
region given by the WMAP data. We use the WMAP parameters for the 
ACDM model (i.e. £^ 0) = 0.73, Q^ = 0.23, fi[ 0) = 0.04, and h = 0.72) 
b . In both models we use the same spectral index n s = 1. The heights 
of the acoustic peaks at small scales (i.e large £) can be affected by the 



b Our data prefers WMAP 3 year data to WMAP 1 year one. 
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following two factors. One is the fact that the scaling of the CDM energy 
density deviates from that of the baryon energy density. Therefore for the 
given CDM and baryon energy densities today, the energy density contrast 
of baryons at decoupling (f^ s ) is getting lower as the coupling is being 
increased. This suppresses the amplitude of compressional (odd number) 
peaks while enhancing rarefaction (even number) peaks. The other is that 
for models normalized by COBE, which approximately fixes the spectrum 
at I ~ 10, the angular amplitude at small scales is suppressed in the coupled 
quintessence. This is shown in the second panel of Figure 1. The third peak 
in this model is smaller than that in the ACDM model. 
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n = 0.001 

n = 0.01 
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Figure 1. (a) CMB large-scale anisotropy power spectra of ACDM (solid line), mini- 
mally coupled n c = (dotted line), and non-minimally coupled n c = 10 — 3 , 10 -2 (dashed, 
dash-dotted line respectively) quintessence models, (b) Same spectra for the entire scales. 



3.2. Matter Power Spectrum 

The coupling of quintessence to the CDM can change the shape of matter 
power spectrum because the location of the turnover corresponds to the 
scale that entered the Hubble radius when the Universe became matter- 
dominated. This shift on the scale of matter and radiation equality is 
indicated 

„(°) 



Pc 



cxp[B c (0o) - B c 



J eq) 



(22) 



where p r and p c are the present values of the energy densities of radiation 
and CDM respectively, and the approximation comes from the fact that 
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the present energy density of CDM is bigger than that of baryons (p c > 
p b ). This is indicated in Figure 2. Increasing the coupling shifts the 
epoch of matter-radiation equality further from the present, thereby moving 
the turnover in the power spectrum to smaller scale. If we define k eq as 
the wavenumber of the mode which enters the horizon at radiation-matter 
equality, then we will obtain 

2tt 



Veq 



(23) 



(i = O.76 

h = 0.70 

Q c = 0.194 

[} b = 0.046 




10 
0.000 



Figure 2. Matter power spectra for the models using the same parameters in Figure 1. 



4. Conclusions 

We have investigated the tracking condition of the quintessence models 
and their stability. We have shown that it is necessary to distinguish the 
tracking condition at the matter dominated epoch and at the radiation 
dominated one. 

We have considered the CMB anisotropy spectrum and the matter power 
spectrum for the non-minimally coupled models. Additional mass and 
source terms in the Boltzmann equations induced by the coupling give the 
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rapid changes of the Newtonian potential $ and enhance the ISW effect in 
the CMB power spectrum. The modification of the evolution of the CDM, 
p c = p c 'a~ 3+ £, changes the energy density contrast of the CDM at early 
epoch. We have adopted the current cosmological parameters measured by 
WMAP within la level. With the COBE normalization and the WMAP 
data we have found the constraint of the coupling n c < 0.01. The locations 
and the heights of the CMB anisotropy peaks have been changed due to 
the coupling. Especially, there is a significant difference for the heights of 
the second and the third peaks among the models. Thus upcoming obser- 
vations continuing to focus on resolving the higher peaks may constrain the 
strength of the coupling. The suppression of the amplitudes of the matter 
power spectra could be lifted by a bias factor. However, a detailed fitting 
is beyond the scope of this paper. The turnover scale of the matter power 
spectrum may be also used to constrain the strength of the coupling n c . 
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